Non-adiabatic scattering of a classical particle in an inhomogeneous magnetic field 
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We study the violation of the adiabaticity of the electron dynamics in a slowly varying magnetic 
field. We formulate and solve exactly a non-adiabatic scattering problem. In particular, we consider 
scattering on a magnetic field inhomogeneity which models scatterers in the composite-fermion 
theory of the half-filled Landau level. The calculated non-adiabatic shift of the guiding center is 
exponentially small and exhibits an oscillatory behavior related to the "self-commensurability" of 
the drifting cyclotron orbit. The analytical results are complemented with a numerical simulation. 



PACS numbers: 73.50.Jt, 03.20.+i, 71.10.Pm 

The dynamics of particles moving in a spatially ran- 
dom magnetic field (RMF) B(r) has been a subject of 
considerable interest in the last few years. The partic- 
ular appeal of the subject is due to its relevance to a 
number of models of strongly interacting disordered elec- 
tron systems in two dimensions (2d). One of the most 
prominent examples is the composite-fermion description 
|jj of a half-filled Landau level. Within this approach, 
the electron liquid in a strong magnetic field is mapped 
- by means of a Chern- Simons gauge transformation - 
to a fermion gas subject to a weak effective magnetic 
field. Precisely at half-filling, the expectation value of 
the Chern-Simons gauge field compensates the effect of 
the external magnetic field. The RMF appears in this 
model after taking static disorder into account: fluctua- 
tions of the local filling factor due to screening of the ran- 
dom scalar potential yield a local mismatch between the 
gauge and external magnetic fields. A number of observa- 
tions ||] of Fermi-surface features near half-filling give a 
strong experimental support to the model of the effective 
magnetic field. Apart from the composite-particle mod- 
els involving fictituous fields, 2d electron systems with a 
real RMF can be directly realized in semiconductor het- 
erostructures by attaching to the latter superconducting 
P| or ferromagnetic [0j5| overlayers. 



The peculiarity of transport properties of 2c? electrons 
in a random field B(r) shows up most distinctly in sys- 
tems with smooth inhomogeneities. The case of long- 
range disorder is most important also experimentally - 
since the compressible state in a half-filled Landau level is 
observed in high-mobility samples. In the latter, a large 
correlation radius of potential fluctuations, d, is fixed by 
a wide "spacer" between the electron gas and the doped 
layer containing ionized impurities. Likewise, inhomo- 
geneities of the RMF created by the ferromagnetic over- 
layers in appear to be fairly long-ranged. In the 
composite-fermion model of the half-filled Landau level, 
the large value of d (as compared to the interelectron dis- 
tance) allows to ignore, in the first instance, the quantum 
interference of scattered waves and describe the electron 
kinetics as purely classical. Quantum localization effects 
become crucial either on exponentially long scales or at 



strong deviations from half-filling, thus motivating us to 
consider a classical particle subject to the RMF. 

This paper is closely related to the previous study |J 
of the electron kinetics in the RMF. The conductivity in 
the field B(r) — B + SB(r) in the limit of strong homo- 
geneous component B was shown in [|| to be completely 
determined by weak non-adiabatic processes in scattering 
on the inhomogeneities 5B(r). A key ingredient of this 
transport theory is the non-adiabatic dynamics on the 
microscopic scale. The purpose of the present paper is 
to study the microscopic processes of the non-adiabatic 
scattering in detail by complementing the analytical ar- 
guments with results of a numerical simulation. 

To this end we formulate a single-scattering problem: 
we introduce a weak homogeneous gradient of the back- 
ground magnetic field and consider the interaction of an 
electron with an "impurity" modeled by a spatially local- 
ized perturbation SB(r). Far away from the impurity the 
electron motion is a slow van Alfven drift along straight 
lines - which are contours of constant magnetic field - 
accompanied by rapid cyclotron gyrations. We consider 
the adiabatic limit, where the shift of the guiding center 
during one cyclotron period, 6, is much smaller than the 
size d of the impurity, d/S 1. In this case, the particle 
continues to drift along the lines of constant B(r), which 
is a manifestation of adiabatic invariance. A weak vio- 
lation of the adiabatic invariance leads to a small shift, 
after the scattering process has taken place, of the mag- 
netic field contour along which the particle continues its 
drift. It is this shift, Ap, that will be studied below, both 
analytically and numerically. The corresponding effect in 
a system with many impurities governs the magnetocon- 
ductivity, as discussed in 

In terms of the complex coordinate z = x + iy, the 
equation of motion at the energy mVp/2 in our scatter- 
ing problem reads z — iQ(z)z, where 
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u> c = eB/mc, B is the background field at y = 0, 
Rc = vf/lu c , the slope e yields a finite velocity of the 
incident particle in x direction, and W(z) is a "scattering 



1 



potential". We assume that e< 1. The guiding center 
coordinate y averaged over the cyclotron orbit, p = (y) c , 
plays the role of an impact parameter. The particle en- 
tering the system at x — — oo with the initial condition 
lmxr^-oo (y) = pi will leave it at x = oo along the tra- 
jectory with (y) c = pi + Ap, where Ap is the desired 
non-adiabatic shift. In order to analyze the scattering 
problem analytically we assume the impurity field to be 
weak, W< 1. Expanding the coordinate z in powers of 
W: z = z + z\ + . . ., we concentrate on the first-order 
term z\. 

The solution y in the absence of the impurity can 
be written in an implicit form for the initial conditions 
z(0) = 0, i(0) = iv F as t(yo) = oj- 1 fl° /R ° dYD~V 2 (Y), 
where D(Y) = 1 - Y 2 (l + \Y) 2 . The coordinate 
yo is confined to the region y_ < y a < y + , where 
y±/R c = (v 7 ! ± 2e— l)/e, and is periodic with the period 
2\t + — t-\ = 2-k/uj, where t± = t(y±). The velocity in x 
direction, xq = —bj c {yQ + ey^/2R c ), may be integrated to 
give xq = (io) c t + £, where £(t) is also periodic with the 
frequency w. The leading terms in the small-e expansion 
are: (x ) c = ^v Fl 

£ = R c [(cosLut - 1) - e(smojt + ^ sm2ut) + 0(e 2 )] , (2) 
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Vo = R c [sinu}t — —e + e(cosujt — — cos2ut) + 0(e 2 )] , (3) 

and lj = uj c + 0(e 2 ). To find yi, we integrate the equation 
of motion for x\ once and use the relation xiio+yiyo = 
following from energy conservation, which yields a first- 
order differential equation for y\ with periodic coefficients 

mm - myi = u c x dt'W[z (t')]y (t') . (4) 

Jo 

For each half-period between the turning points t = t±, 
at which yo changes sign, the solution of Eq. (4) can be 
represented in the form 

yi(t)=w c yo(t) tdt'K(t') [ dt"W[z (t")]y (t") . (5) 
Jo Jo 

Here the factor K(t) = i (i)/j/g(t) behaves singularly as 
(t — t±)~ 2 in the vicinity of t±. Matching two branches of 
yi(t) across the turning point requires that the contour 
of the t' integration in Eq. (5) be displaced into the com- 
plex plane so as to pass round the singularity on the real 
axis [since the term (f — 1±)~ 1 is absent in the integrand, 
the contour can be shifted in either half- plane] . This cor- 
responds to the contraction f t *^ dt(t — t±)~ 2 — ► — 2t _1 
at t — > 0. With this choice of the contour of integration, 
Eq. (5) gives the sought solution at all t. 

To extract the non-adiabatic shift from Eq. (5), we 
observe that the integral over t" converges to a con- 
stant at large t' while (yo) c = 0, so that the average 
= (yi(t -» oo)) - (yi(t -> -oo)) is given by 
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where a(e) — u c (yo{t) § dt 1 K(t')j is expressed in 

terms of the unperturbed solution. At e — > the con- 
stant a — > — 1. Note that, apart from the shift Ap, the 
asymptotics of Eq. (5) at large t contains an oscillating 
term — eluj c t cos u) c t, whose amplitude diverges linearly 
with growing t. This divergency is an artefact of the per- 
turbation expansion in W and reflects the fact that the 
shift Ap is accompanied by the change of the frequency 
(eto c / R c )Ap. This asymptotics could have been equiva- 
lently used to find Ap = —I. 

In evaluating / we first assume, for simplicity, that 
W{z) depends on x only, which allows to expand I as 

dtJ2^ d 2 w &Ft)C(t)y (t). (7) 

' n 

Since W(^vpt) is a smooth function of t on the scale 
of uj^ 1 , the leading contribution to I comes from taking 
the first harmonic of the integrand - higher harmonics 
will involve exponentially smaller Fourier components of 
W(^vpt). This yields, after integration by parts and re- 
summation, 
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The integral in Eq. (9) can be evaluated at e 1 by the 
saddle-point method to give 



Ap = —2vp\ — cos 

V 7T V e 4 



(10) 



J dt COS 0J c t W {^°Ft — Rcj 



This equation is a paramctrically exact solution of the 
scattering problem at e — ► 0. It expresses the non- 
adiabatic shift in terms of the asymptotics of the Fourier 
transform of the smooth function W - thus demonstrat- 
ing explicitly the exponential smallness of Ap. The pa- 
rameter which governs the exponential falloff of Ap is 
d/5 ^> 1, where d is a characteristic scale of variation of 
W and S = TreR c is the shift of the guiding center af- 
ter one period of the cyclotron rotation, while the ratio 
d/R c may be arbitrary. Though the limit of a smooth 
inhomogencity on the scale of R c is historically most 
closely associated with the notion of adiabaticity, the 
parameter d/R c plays no role in Eq. (10). It is worth 
noting that the pre-exponential factor given by the first 
line of Eq. (10) is a non-analytic function of the small 
parameter e and, moreover, happens to oscillate wildly 
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as e — > 0. These oscillations are due to the commen- 
surability of two length scales R c and S. Remarkably, 
the series of the geometric resonances which constitute 
the oscillations is associated with the properties of the 
unperturbed solution ( "self-commensurability" ) and not 
with the shape of the scatterer. Note that to get the os- 
cillations one has to sum up all terms in the expansion 
(7) even if W is smooth on the scale of R c . Another 
peculiar feature of the non-adiabatic shift is related to 
its sensitivity to the phase 4> of the cyclotron rotation 
of the incident electron. Specifically, changing the po- 
sition of the scatterer by Ax leads to the oscillations of 
Ap(<j>) = Ap m cos(0-0 o ) oc dtco&(uj c t+(j))W(^v F t), 
where (p = 2ttAx/5, so that the shift vanishes periodically 
with varying initial conditions. 

In the composite-fermion problem, a charged impurity 
located the distance d from the plane occupied by the 
electron gas creates the axially symmetric perturbation 
W(y) = VF d 3 [(r-R) 2 +d 2 ]~ 3 / 2 . We choose the impurity 
position R = (0,pi), so that pi has the meaning of the 
impact parameter with which the guiding center is inci- 
dent on the impurity. Extracting the first-order Fourier 
component of W[zo(t)] in the same way as in Eqs. (8), (9), 
we get 
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The function W in this equation has branch points at 
H = t' s (t), where t' s (t) = ±^#+W*FftF, which 
determines the exponentially small value of the integral 
over t'. Specifically, the second line in Eq. (11) reads 
V^W (2nd/S) 3 (^\t' s (t)\)- 3 / 2 cxv(-uj c \t' s (t)\). The re- 
maining integral over t can be done by the saddle-point 
method. The cumbersome general expression reduces to 
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in the limit of a long-range impurity potential, when 
d, \pi\ > R c . Here d = (d 2 + p 2 ) 1/2 . In the opposite 
case, when R c is larger than both d and \pi\, we get 



d 2 2 ( d 

Ap ~ 8tt Wq — r cos 4> cos - exp — 2ir— 
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(13) 



Now let us turn to a numerical simulation. A typi- 
cal trajectory resulting from the numerical integration of 
Eq. (3) is shown in Fig. 1. Though the pertubation of 
the drift trajectory is seen to be large in this particular 
example, the non-adibatic shift is almost invisible. The 
inset to Fig. 2 shows the periodic dependence of the shift 
on the phase <f>, which agrees with Eq. (12). The main 
panel clearly demonstrates the adiabatic character of the 
scattering at d/S ^> 1: the magnitude of the oscillations 
Ap m plotted against S^ 1 is seen to fall off exponentially, 
also as predicted by Eq. (12). The slope c of the exponen- 
tial decay InAp = —2ircd/5 is found to be 1.3, which is 
somewhat larger than the value c = 1 following from Eq. 
(12). This is because the perturbation of the drift tra- 
jectory is not negligible in this case. No oscillations with 
changing e could be reliably seen in Fig. 2, which should 
also be ascribed to their smearing due to the curvature of 
the trajectory. Indeed, the commensurability of R c and 5 
cannot be maintained in the whole region of interaction 
if the drift trajectory is strongly perturbed, as clearly il- 
lustrated by Fig. 1. The non-adiabatic shift at a smaller 
amplitude of the interaction Wo is shown in Fig. 3. In 
this case, the perturbation of the trajectory was weak, 
and the particle was essentially drifting right "through 
the impurity" with a small impact parameter pi < d, R c . 
Decreasing Wo narrows the range of the numerical simu- 
lation but clearly reveals the oscillatory behavior of Ap m 
with changing e. The exponential dependence in Fig. 3 is 
fitted very well by e r 2lxd l & (i n accordance with Eq. (13)) 
with no fitting parameter in the exponent. The oscilla- 
tions of Ap m with changing R c /S (inset of Fig. 3) have a 
period close to 1/2 and are in good agreement with the 
analytical prediction Ap m oc | cos 27r ^ c | . 




FIG. 1. Typical trajectory of a particle scattered on a 
"magnetic impurity" [Eq. (1)]. The lines of constant magnetic 
field -B(r) are shown. The strength of the impurity Wo = 2.2. 



The last three factors in Eqs. (12), (13) reflect the features 
of the non-adiabatic shift discussed above: the exponen- 
tial smallness, the oscillations with changing e, and the 
oscillatory dependence on the phase (p. Note that all the 
formulas above imply that the drift trajectory is only 
slightly perturbed by W(r). 
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FIG. 2. Amplitude of the non-adiabatic shift Ap m as a 
function of S = 2nVd/u} for a strong impurity: Wo ~ 2.2. The 
solid line is a fit Ap m oc exp[— 1.3 x 2ird/S\. Inset: Shift Ap as 
a function of the phase <f) for different values of the parameter 
2Tvd/6 = 4.8(0), 5.8(D), 7.8(0). 



which models a scatterer in the composite-fermion theory 
of the half-filled Landau level. The non-adiabatic shift of 
the guiding center is exponentially small and exhibits os- 
cillations with i? c /<5, where R c is the cyclotron radius, S 
the shift of the guiding center along the drift trajectory 
after one cyclotron period. The oscillations are related 
to the "self-commensurability" of the drifting cyclotron 
orbit. The analytical results are in full agreement with 
the numerical simulation. A detailed numerical study 
of the transport properties of the fcrmions in a random 
magnetic field will be presented elsewhere. 
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FIG. 3. Non-adiabatic shift for a weak impurity 
(Wo = 0.07) at pi/d = 0.4 - 0.9 and large R c /d = 1.5 - 2.5. 
The solid line is a fit Ap m oc exp(— 27rd/<5). Inset: Oscillations 
of Ap m as a function of R c /S. 

In conclusion, we have studied the violation of the 
adiabaticity of the electron dynamics in a slowly vary- 
ing magnetic field. We formulated a scattering problem 
which has been solved exactly. As a particular example 
we considered scattering on a single "magnetic impurity" 
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